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Abstract. Active dumbbell suspensions constitute one of the simplest model sys- 
tem for collective swimming at low Reynolds number. Generalizing recent work, we 
derive and analyze stroke-averaged equations of motion that capture the effective 
hydrodynamic far-field interaction between two oscillating, asymmetric dumbbells 
in three space dimensions. Time- averaged equations of motion, as those presented 
in this paper, not only yield a considerable speed-up in numerical simulations, 
Q ' they may also serve as a starting point when deriving continuum equations for 

C/5 , the macroscopic dynamics of multi-swimmer suspensions. The specific model dis- 

cussed here appears to be particularly useful in this context, since it allows one 
2 ' to investigate how the collective macroscopic behavior is affected by changes in 

the microscopic symmetry of individual swimmers. 

^ ! 1 Introduction 

The rich collective behavior of active biological systems, such as flocks of birds, schools of fish 
or bacterial suspensions, has attracted considerable interest in biophysics in recent years [TJ 
[2l[3l|4j[5l[6] . Due to the complex multi-scale nature of the observable patterns [7j|8] , theoretical 
approaches are manifold, ranging from discrete 'microscopic' descriptions, that account for the 
dynamics of individual organisms [2,4,5J, to field-theoretic continuums models, that aim to 
^■f-N , capture dynamical features and symmetries on the macro-scale [9l ll0l l6]. The arguably most 

f^ ' important problem in this context is to understand which microscopic properties determine (or 

are reflected in) the macroscopic behavior [U|B]. 
i— V . Generally, linking micro- and macro-behavior is mathematically difficult as it typically in- 

f— ^ ' volves both spatial and/temporal averaging (coarse-graining) of nonlinear equations of mo- 

tion [TT1IT21IT3] . Valuable insight can be gained by considering simplified model systems that 
allow for explicit testing of averaging procedures |14[|15| . In the present paper, we intend to 
study the validity of temporal coarse-graining for a simple model |111I14) of collective swim- 
ming at low Reynolds number, a regime relevant to the motion of bacteria and small algae |16l 
;h ' [T71IT5] . Specifically, generalizing recent work [T^ITTIIH] . we derive time-averaged equations of 

5t 1 motion that govern the effective, three-dimensional hydrodynamic interaction between actively 

oscillating, asymmetric dumbbelQ pairs. 

It has been known for a long time [16] that - owing to the time-reversibility of the Stokes- 
equation - isolated, force-free dumbbells cannot swim in the zero- Reynolds number limit. How- 
ever, when the phases of two identical, periodically oscillating dumbbells at finite separation are 
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detuned, then collective swimming by m.eans of hydrodynamic interaction becomes possible |11) . 
Thus, dumbbell 'swimmers' may be regarded as the simplest model for collective swimming at 
low Reynolds number. Related recent work focusses either on the one-dimensional case [11^ or 
on symmetric two-sphere swimmers [H]. Here, we generalize the analysis of Ref. [13] to the 
asymmetric case. The motivation for this is as follows: 

Symmetric dumbbells do not possess an intrinsic orientation. Consequently, a corresponding 
field-theoretic continuum description must be based on the Q-tensoo since their local mean 
orientation (polarization) field is trivially zero, 11 = 0. By contrast, asymmetric dumbbells 
may possess a nontrivial mean orientation field, n(i,x) ^ 0, which in d > f space dimensions 
can change due to hydrodynamic interactions. Hence, asymmetric dumbbells appear to be the 
simplest model system for studying how the symmetry of the microscopic constituents affects 
both hydrodynamics and orientational order in the continuum limit. 

However, before one can derive the corresponding field equations by means of standard 
methods [5] , one first has to identify the effective equations of motion for a dumbbell's position 
and orientation change by averaging the hydrodynamic interactions with the other swimmers 
over a swimming stroke. It is, therefore, the purpose of the present paper to (i) provide explicit 
expressions for the time-averaged, effective interaction forces, and to (ii) verify their validity 
by comparing the exact microscopic with the coarse-grained dynamics. 



2 Stroke-averaged equations of motion 

We first summarize the microscopic equations of motions of the dumbbell model. Subsequently, 
the corresponding stroke-averaged equations will be discussed. 



2.1 Microscopic model 

We consider a system of 5 identical dumbbells. Each dumbbell swimmer consists of two spheres, 
with radii a\ and a2- At low Reynolds numbers, inertia is negligible and the state of the system 
at time t is completely described by the spheres' position coordinates {Xq} = \Xi^^\(t)\ with 
a — I, . . . , 25* labeling the spheres, and « = 1, 2, 3 the space dimension (throughout, we adopt 
the Einstein summation convention for repeated Latin indices). Neglecting rotations of the 
spheres, their dynamics is governed by the overdamped equations [TTlfTniin] 

25 

^(m)(i) = 2J '^{oti){fi2)F(fij) (1) 

,3 = 1 

where X := dX/dt is the velocitjo The hydrodynamic interaction tensor "K couples the de- 
terministic force components i^(/3i) that act on the individual spheres. Generally, the vector 
F — {^(^j)} may comprise contributions from internal forces, i.e., those required to bind two 
spheres to form a dumbbell, as well as from external force fields (gravity, etc.); however, in the 
present paper, we shall assume that external forces are negligible. 

In our numerical simulations, "K is given by the Rotne-Prager-Yamakawa-Mazur tensor [THl 
[2nl[2T1[22l[23] 

^^'^'-''^'^^^ - - - ^^^ (2) 



Q-tensor := second moment tensor of orientations [^ 

Ref. [14] discusses how to include thermal fluctuations in Eq. fl}. 
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where r^^i := Xai — xpi, a ^ (3, and rap '■— |xq — x^|. Analytical formulas presented below 
are based on an Oseen approximation, which neglects the r~g term in Eq. ([3]). The diagonal 
components ^ describe Stokesian friction in a fluid of viscosity fi. The off-diagonal compo- 
nents ([3]) model hydrodynamic interactions between different spheres. Note that 5{ is positive 
definite for rap > Oq. -I- a^ and divergence-free, J^a ^{I3j)^iai)(i3j) = with 9(^j) :— d/dx^jj^y 

We still need to specify the intra-dumbbell force F. Let us consider the dumbbell a, formed 
by spheres a — 2cr — 1 and (3 = 2cr, and denote its length by d°'{t) := |X^(t) — XQ,(t)|. Neglecting 
external force fields, we shall assume that the two spheres are connected by a harmonic spring 
of variable length 

L''{t)=e + Xsin{ujt + ip''), £> aa + ap + X. (4) 

In this case, -F(^i) = —d(pi)U where 

C/ = ^f/^ C/-(i,d-) = ^K-L-(ir. (5) 

The dumbbell swimmer is called passive if the stroke amplitude is zero, A = 0, and active if 
|A| > 0. As discussed below, the phase parameter (/3°' is important for the interaction between 
two or more dumbbells. 

For the overdamped description ([T} to remain valid, the driving must be sufficiently slow. 
More precisely, we have to impose that T^ ^ Tq ^ T^, where T^^ := I-kJuj is the driving 
period, Tq '■— 27r/y^fco/Ma the oscillator period for a sphere of mass Af^, and T^ :— Ma/ "fa 
the characteristic damping time. This restriction ensures that the dumbbells behave similar to 
shape-driven swimmers, i.e., df ~ L'^{i) is a useful approximation in analytical calculations. 



2.2 Coarse-grained mesoscopic dynamics 

We next summarize the stroke-averaged equations of motions for the dumbbell positions and 
orientations, obtained by applying the procedure described in the Appendix of Ref. 14 . Below, 
the resulting effective equations of motion will be compared with numerical simulations of the 
microscopic model equations ((!]). 

Each dumbbell can be characterized by its orientation vector 

N-W^ ^^-'^^^Y ^ = 1,---,^ (6) 

and a suitable position coordinate 

R"(0=/32X2.+/3iX2._i, /3i+/32 = l, /3i/2>0. (7) 

For example, the choice Pi = 132 — 1/2 corresponds to the geometric center |llj 

R^(t):=i(X2.+X2.-i), (8) 

Here, we shall consider j3i — ai/{ai + 02), defining the center of hydrodynamic stress [22l [6| 

RS(t) :- «^X2.+aiX2.-i ^ ^g^ 

ai + 02 

For very small asymmetries ai ~ 02, the geometric center Rq practically coincides with the 
hydrodynamic center R|^. For strongly asymmetric dumbbells with ai <^ 02 or oi 3> 02 the 
hydrodynamic center R^ presents the more appropriate choice, as it is the 'slower' variable. 
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The basic idea of the stroke-averaging procedure [TTlfT^fT^] is to focus on the dynamics of 
averaged position and orientation coordinates R(i) and N'^(t), defined by 



J. I^rxi fry -f _1_T^ /O 

N-(t) := i / ds N-(s) , R-(i) := i / ds R-(s). (10) 



Here T = 2-k juj denotes the period of a swimming stroke. If N'^(i) and R'^(i) are slowly varying 
functions of time, one can approximate 



N'^-N , R'^-R, -/ ds/(N-(s),R-(s))~/(N-(t),R-(t)) (11) 

^ Jt-T/2 



for any sufficiently well-behaved function /. For nearly symmetric dumbbells, the approxima- 
tions (jlip are justified for both Rq and Rh, whereas for strongly asymmetric ones they usually 
only hold for the hydrodynamic center Rjj. 

Using Eqs. (|lip. one can derive from the microscopic model equations H]) the correspond- 
ing deterministic stroke-averaged equations of motions pilfT^fTB] . by assuming that: (i) The 
dumbbells are force-freeQ and approximately shape driven, i.e., d" := |X2cr — X2cr-i| — L'^{t)- 
(ii) The dumbbells are slender, i.e., sphere radii ai/2 and stroke amplitude A have about the 
same size, but are much smaller than the dumbbell's mean length t. (iii) The ensemble is dilute, 
meaning that the distance D"^,^ := jD^'^jjl := IRq/jj — Rq/hI between dumbbells a and p is 
much larger than i. 

Adopting the simplifications (i)-(iii) and restricting to two-body interactions, one finds the 
following coarse-grained equations of motion for the hydrodynamic center, 

^- = fiTl)Er+E^r, (12) 




N! = ^{5., - ivrK) { ( ^f^ ) ^ K' + E ^r ^ (13) 

PT^cr P#0" 

where the stroke-averaged hydrodynamic interaction terms to leading order in \j t are given by 

^.— <.-.".S(^)'(^)(4i)^ 

iVf(l-3r2-3s2-6gsr+15s2r2), (14) 

jr^^^...(.'-^')|;Q)'(jj^)'x 

{iVf (2s + 4gr - lOsr^) + i)^^(l + 2^^ - Ss^ - Sr^ - 20gsr -|- SSsV^)}, (15) 



Here, 



-r-™(.-.'.|Q)^(f)(^)^ 

L>Hfc(l + V - 5s2 - 5r2 - 20gsr -I- SSsV^), (16) 

^Hfe (3s + 6rg -I- Gsg^ - Ts^ - 2\sr^ - 42qs^r + BSsV^) . (17) 

yl :=: 20102/(01+02) (18) 



* If the internal forces required contract the dumbbells are central forces, then the force-constraint 
implies that the torque-free constraint is automatically fulfilled. 
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denotes the harmonic mean of the sphere radii, the unit vector D^'' ;= T)'^ /\'D'^\ gives the 
orientation of the distance vector D^'' — R^ — R^, and s, r, q abbreviate the projections 

s:=DS^7V;, r:=D^PN^, 9 :- A^/A^/- (19) 

One readily observes two prominent features: Firstly, the stroke-averaged interactions terms 
/, J, K, L vanish if the phases ip'^ and (fP differ by multiples of tt. Secondly, in the dilute limit 
|D^''| 3> i the leading contribution to R^ is given by the /-terms, which decay as \T)'^\~^ [TI] . 
By contrast, for symmetric dumbbells with ai = 02 the first term in Eq. (IT^ is absent and one 
recovers the equations recently derived in Ref. |12) : hence, in this case, the effective interaction 
force decays asymptotically as |D^''|~*. Similarly, asymmetric dumbbells experience changes in 
rotation N"' which decay as |D^''|^'^, whereas the effective rotation interaction for symmetric 
dumbbells decays as jD^''!"^. 



3 Numerical tests 

We verify the validity of the stroke-averaged equations of motion (TT^ by comparing with time- 
resolved simulations of the microscopic dynamics ([T|). To this end, we numerically simulate 
both microscopic and coarse-grained equations of motions using a CUDA algorithm that is 
described in detail in Ref. [13]. We first compare our results to those of Lauga and Bartolo [IT] 
for the one-dimensional case. Our findings agree with theirs for most setups, they differ for one 
particular configuration (see discussion below). We conclude our numerical analysis by testing 
the coarse-grained equations of motions for various three-dimensional dumbbell configurations. 
Generally, when comparing symmetric and asymmetric dumbbell configurations, one should 
keep the quantity A — 2aia2 / {ai + a2) constant, in order to keep the 'symmetric' J-contribution 
in Eq. ^ fixed. 



3.1 One-dimensional case 

We first consider aligned dumbbell pairs (see Fig. [IJ as studied by Lauga and Bartolo 1111 . 
Aligned dumbbells do not change their orientation and Eqs. p^ reduce to 



A 



2 



3 

' '.crp 



where R^ denotes the coordinates along the common axis. The lines in Figs. [T] (a) and (b) 
represent the dynamics of aligned dumbbell pairs as predicted by Eq. (1201) , and symbols indicate 
the results of corresponding microscopic model simulation. Following Lauga and Bartolo [TT], we 
quantify collective motion of the dumbbell pairs in terms of their mean collective displacement 
(solid lines/filled symbols in Fig[T]), 

Wit) = liRUt) + Rhit)l (21) 

and their mean relative distance (dashed lines/unfilled symbols in Fig [T]) , 

ARll{t)^RUt)-Rhit)- (22) 

The quantity R^(t) characterizes the net motion of the dumbbell pair, whereas AR^{t) indi- 
cates whether the dumbbells move towards or away from each other. 
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Fig. 1. Depending on the symmetry of their initial configuration, aligned dumbbell pairs ex- 
hibit qualitatively different collective motions. In both diagrams, lines were obtained by nu- 
merical integration of the stroke-averaged equations ()20|) . whereas symbols show the simula- 
tion results for the microscopic spring-based dumbbell model ((T]). Solid lines and filled sym- 
bols depict the mean displacement '^{t) -^{Q) = {[R\(t) + R\i{t)] -[R^{Q) + R\i{Q)]}/2 of 
the hydrodynamic centres. Dashed lines and unfilled symbols indicate the relative separation 
AR^{t) - Z\i?H (0) = [i?H(f) - -Rh(*)] - [-Rh(O) - R\i{Q)]- (a) For mirror symmetric configurations, the 
mean distance between the dumbbells remains constant and the dumbbell pair can move in either di- 
rection. In particular, asymmetry can enhance the collective speed, see triangles, (b) For translation 
invariant configurations, asymmetry does not affect the collective pair velocity, but depending on the 
initial orientation the dumbbells either approach each other (red triangles) or move away from each 
other (black diamonds). Simulation parameters are comparable to those of Lauga and Bartolo [11| : 
Initial separation AR^{Q) — -Rh(O) — ^h(O) ~ 10^, mean dumbbell length t = 5/.im, driving fre- 
quency u) — SOOs"^ (time is given in units of the stroke period T = 2-k /u), stroke amplitude A = Q.ll , 
A = 2aia2/(ai + 02) — 0.15£ with 02 — ai/2 for asymmetric dumbbells, phase difference tp^ —ip^ = 7r/2. 
For the microscopic model: spring constants fco — O.OOlkg/s^, viscosity ^ — IQ-^ kg/(ms), particle mass 
density g = 10'^ kg/m''; simulation time step At « IQ-^T. We note that the quantitative difference 
between symbols and lines decreases when choosing smaller ratios (aa/l), cf. discussion in Sec. 4.1 of 
Ref. [H]. 



As evident from the two diagrams in Fig. [U we find qualitatively different behavior depend- 
ing on the symmetry of the initial dumbbell pair configuration. Figure [1] (a) summarizes results 
for reflection symmetric initial conditions. In this case, the pair can move in either direction, 
depending on whether the smaller spheres point towards or away from each other, but the mean 
distance between the two dumbbells remains constant. By contrast, the translation invariant 
configurations in Fig. [1] (b) always move in the same direction with nearly identical speeds 
for symmetric and asymmetric dumbbell pairs. We note that our results for the translation 
invariant setups essentially agree with those obtained by Lauga and Bartolo [TT]. However, we 
find a different behavior for mirror symmetric configurations, cf. Fig. 2 (a) of their paper [llj . 
In particular, our results show that under suitable conditions collective motions can be signif- 
icantly enhanced by asymmetry, see red triangles in Fig[T] (a). The fact that asymmetry can 
increase effective hydrodynamic interactions could be of relevance for the construction of more 
efficient micropumps [24j . 

Furthermore, Fig. [3] shows how both the collective displacement and the relative separation 
per period vary with the distance between the dumbbells. As in Fig. [TJ lines indicate the 
prediction based on the stroke-averaged equations (J12p . whereas symbols indicate the simulation 
results for the related spring-based model (symbols, colors, and line-styles refer to the same 
configurations/model parameters as in Fig. [U respectively) . Remarkably, the dynamics of the 
microscopic model is very well described by the averaged equations (|12p down to distances 
of a few body lengths. This is a bit surprising given that the stroke-averaged equations are 
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Fig. 2. Distance dependance of the collective motion and separation for aligned dumbbell pairs. Line 
styles and symbols correspond to the same configurations and simulation parameters as used in Fig. [T] 
For a refiection symmetric configuration such that the small spheres point outwards, see black diamonds 
in diagram (a), the collective velocity reverses its sign at a distance of approximately 11^. Qualitatively, 
this change is correctly refiected by the stroke-averaged equation (|12|1 . see solid black line in diagram 
(a). Remarkably, the stroke-averaged equations describe the microscopic dumbbell dynamics reasonably 
well down to distances of a few body lengths. 



based on a far-field expansion. In particular, as also correctly predicted by the stroke-averaged 
equations (J12p . for the mirror symmetric configuration with small spheres pointing outwards 
the pair velocity reverses its sign at a distance of approximately 11£; see black diamonds in 
Fig.Ka). 



3.2 Three-dimensional case 



In two- and three-dimensional systems, hydrodynamic interactions can not only lead to collec- 
tive translational but may also induce orientational changes. To illustrate both effects and to 
verify the validity of the angular parts in Eqs. (I12p - (I17I) . we consider various different initial 
configurations as sketched next to the diagrams in Fig. [3] and 01 In each of the diagrams, blue 
symbols/squares indicate the change of the relative orientation of the two dumbbells, quantified 
through the change of the projection 



Aq{T) := g(0) - q{T) , 



g(t):=7Vf(t)7Vf(i), 



(23) 



where N°'(i) and N''(i) denote the orientation of the dumbbells a and p at time t, respectively. 
It is important to note that, for symmetric dumbbells with a\ = 02, the orientation vector 
is not uniquely defined (in this case, the coarse-grained equations of motions are invariant 
under the transformation N 1-^ ~N). For asymmetric dumbbells, however, the orientation 
can be uniquely characterized by means of the different sphere radii. In our plots, we fix the 
orientation as pointing from the filled to the unfilled sphere in all cases. 

Furthermore, the red lines/triangles show the change in the relative separation of the two 
dumbbells after one period. 



AD^{T) := [DuiO) - Dh(T)] /£ , Dait) := |R5(t) - R£(t)|/1 



(24) 



For AD}i{T) > the dumbbells approach each other; for AD^iiT) < they move away from 
each other. 

To illustrate the effects of a change in symmetry, we plot the results for symmetric and asym- 
metric dumbbell pairs next to each other. As one may readily observe, asymmetry does strongly 
affect both the translational and orientational motions of the dumbbells. This suggests that ge- 
ometry can play an important role for the emergence of orientational (dis)order in suspensions 
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Fig. 3. In-plane rotation of the initial swimmer positions. Symbols are obtained by numerical 
simulations of the microscopic equations of motions and lines indicate the results of the stroke- 
averaged dynamics. Blue solid lines/blue squares denote the change of the relative orientation Aq 
over a stroke period. Red dotted lines/red triangles indicate the change of the relative distance ADh 
over a stroke period. Swimmer parameters are identical to those used in the Id tests except for 
A = 2aia2/{ai + a2) = O.U. 

of hydrodynamically interacting organisms. Generally, Figs. [3] and S] confirm the coarse-grained 
equations of motions (|12I) - P7)) are able to quantitatively reproduce the characteristic features 
of the microscopic dynamics. 



4 Summary 

Active dumbbell suspensions present a useful model for studying collective motion and orien- 
tational ordering due to hydrodynamic interactions at low Reynolds number. In this paper, 
we have discussed the coarse-grained equations of motion for the effective hydrodynamic in- 
teraction between actively oscillating, asymmetric dumbbell pairs. Our analysis shows that 
the stroke- averaged equations p^ - (IT7)) are able to capture the main features of the time- 
resolved, microscopic dynamics at moderate-to- low densities (intermediate-to-large distances). 
Thus, equations of the type (|12p - (jl7p provide a convenient mesoscopic description which, for 
example, can be used as a starting point for the derivation of coarse-grained macroscopic field 
theories [lO] . 

The coarse-grained equations of motions discussed above can be readily implemented in 
GPU-based CUDA simulations analogous to those described in Ref. [H]. The good agreement 
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Fig. 4. Out-of-plane rotation of the initial orientation of the second swimmer by an angle (/f>a around the 
indicated axis. Symbols are obtained by numerical simulations of the microscopic equations of motions 
and lines indicate the results of the stroke-averaged dynamics. Blue solid lines/blue squares denote the 
change of the relative orientation Aq over a stroke period. Red dotted lines/red triangles indicate the 
change of the relative distance ADn over a stroke period. Swimmer parameters are identical to those 
used in the Id tests except for A = 2aia2/(ffli + 02) = Q-U. 

between the averaged dynamics and the microscopic model simulation confirms that our CUDA 
algorithm }14] works correctly even at relatively low densities, when hydrodynamics interac- 
tions effects are relatively weak and algorithms may become prone to numerical instabilities. 
With regard to the future, we hope that the combination of GPU-based simulation techniques 
and systematic coarse-graining will enable us to quantitatively compare particle simulations 
with continuum field theories. The asymmetric dumbbell model considered here appears to be 
particularly promising in this context as it allows one to investigate how microscopic symmetry 
breaking affects macroscopic behavior. 
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